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O ■ Abstract 

\o 

• We solved the Schrodinger equation for a particle in a uniform magnetic field 

in the n-dimensional torus. We obtained a complete set of solutions for a 
» ^ ■ broad class of problems; the torus T n = R n / A is defined as a quotient of 

■ the Euclidean space R n by an arbitrary n-dimensional lattice A. The lat- 

tice is not necessary either cubic or rectangular. The magnetic field is also 
arbitrary. However, we restrict ourselves within potential-free problems; the 
J^~j ■ Schrodinger operator is assumed to be the Laplace operator defined with the 

j_j ■ covariant derivative. We defined an algebra that characterizes the symmetry 

of the Laplacian and named it the magnetic algebra. We proved that the 
space of functions on which the Laplacian acts is an irreducible representation 
space of the magnetic algebra. In this sense the magnetic algebra completely 
characterizes the quantum mechanics in the magnetic torus. We developed a 
new method for Fourier analysis for the magnetic torus and used it to solve 
the eigenvalue problem of the Laplacian. All the eigenfunctions are given in 
explicit forms. 
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1 Introduction 



In this paper we solve the Schrodinger equation for a particle in a uniform magnetic field in 
an n-dimensional torus. The problem looks plain at first sight but actually it turns out to 
be a hard problem. Hence we begin this paper by a quick explanation of the problem. After 
that we will describe our strategy to solve it. Subsequently we will briefly mention studies 
by other people and describe our motivation of this study. At the end of Introduction we 
will give guides for quick access to main results of this paper. 

An n-dimensional torus, or n-torus, is defined as T n = R n /Z n . In the coordinate a 
point (t 1 , . . . , V + 1, . . . , t n ) is identified with (t 1 , . . . , P , . . . , t n ) for each j = 1, . . . , n. The 
eigenvalue problem of the ordinary Laplacian in the torus is the equation 

-a/ = -e(^) /-/ (1-1) 

with the periodic boundary condition 

/(/'' /•' • i n /(/' /•' /")• (i.2) 

The eigenvalue problem can be immediately solved by Fourier expansion. A plane-wave 
function 

Xk (t\...,t n ) = e 2m ^U k ^ (1.3) 

with quantized momenta kj G Z is a solution. The whole set of eigenfunctions {xk \ (ki, ■ ■ ■ , 
k n ) G Z n } constitutes a complete orthonormal set of the space of periodic functions over 
the torus. This is a well-known result. 

In this paper we would like to solve an eigenvalue problem of the magnetic Laplacian; the 
magnetic Laplacian is defined by replacing the partial derivative in the ordinary Laplacian 
by a covariant derivative as 

A/= £^-2^) (A- 2,a,)/. (1.4) 
Here Ai is a component of the U(l) gauge field 

^ = ^E^-|* , ' + ai (1-5) 

with integers {(f) ji = —<pij} and real numbers {otj}- The gauge field A = J2?=i M dt l generates 
a uniform magnetic field B = dA = (l/2)J2 1 ji=i4>jidP A dt l . Moreover, we would like to 
consider a general oblique torus T n = R n / A; A is an n-dimensional lattice. Edges of the unit 
cell of the lattice do not necessarily cross at a right angle and they do not necessarily have 
a same length. Hence we introduce a metric g^ 1 in the definition of the magnetic Laplacian 
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(11-4)1 to take inclined and stretched or shortened unit cells into account. The eigenvalue 
problem of (jl.4j) is accompanied by the condition 

f(t\ ..../•''• 1 /") e^^Li *»*f(t\ t n ), (1.6) 

which we call a twisted periodic condition. Thus the plain problem (|1.1|) with (|1.2|) is gener- 
alized to the 'magnetic' problem (jl.4|) with (|1.6j) . At first glance it looks rather straightfor- 
ward to generalize the problem in this way but it is actually highly nontrivial and difficult 
to generalize the solution. 

Let us see where the difficulty lies. In the case of the ordinary Laplacian, the plane-wave 
solution (jl.3)) is a simultaneous eigenfunction of the momentum operators 

d 

Pj = - 1 -^ (j = l,...,n) (1.7) 

as PjXk = 2'xkjXk- The Laplacian can be expressed in terms of the momentum operators as 
—A = Y^j=i{Pj) 2 an d, of course, it commutes with the momentum operators. Thus integers 
(fci,..., k n ) are 'good' quantum numbers. Then the whole set of simultaneous eigenf unctions 
{Xk} forms the complete solutions of the Laplacian problem. This is the way how Fourier 
analysis works. However, when we turn to the magnetic Laplacian, we may seek for a 
simultaneous eigenfunction of 'magnetic' momentum operators 

Pj = -i(j--2-KiA^J (j = l,...,n). (1.8) 

But such a simultaneous eigenfunction does not exist because magnetic momenta do not 
commute with each other and instead exhibit commutators 

[pj,Pi] = 2iri</>j l . (1.9) 

The magnetic Laplacian can be still expressed in terms of the magnetic momentum operators 
as —A = J2j=i{Pj) 2 but it does not commute with pj. Hence, the strategy of ordinary Fourier 
analysis does not work well for the magnetic Laplacian. 

To solve the problem we developed a new method, which we call Fourier analysis for 
the magnetic torus. This is a main subject of this paper. Let us describe our strategy: 
First, we will find a group of operators that commute with magnetic momentum operators. 
We call the group a magnetic translation group. Second, we enlarge a family of operators 
to define an algebra, which includes magnetic momenta and magnetic translations as its 
elements. We call the algebra a magnetic algebra and construct its representations. Third, 
we show that the space of twisted periodic functions over the torus is actually an irreducible 
representation space of the magnetic algebra. By diagonalizing a maximal commutative 
subalgebra of the magnetic algebra we obtain a complete orthonormal set of twisted periodic 
functions. This set of orthogonal functions provides a kind of unitary transformation as the 
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set of plane-wave functions provides the Fourier transformation which bridges between the 
momentum space and the real space. We note that it is easy to diagonalize the Laplacian 
in the momentum space. Finally, we get a whole set of eigenfunctions in the real space by 
applying the unitary transformation. In this procedure the third step is the hardest part 
and is actually accomplished by lengthy cumbersome calculations. However, the strategy is 
clear. 

We would like to briefly review studies by other people on spectral analysis in magnetic 
field. Brown Q first examined the symmetry structure of the Schrodinger equation for an 
electron in a lattice in a uniform magnetic field and found that the symmetry is described by 
a noncommutative discrete translation group. At the almost same time Zak [SJ also found 
the same symmetry structure and named the group a magnetic translation group (MTG). 
Zak (2] immediately built a representation theory of the MTG in the three-dimensional 
lattice. From the viewpoint of functional analysis, Avron, Herbst, and Simon have been 
studying spectral problems of the Schrodinger operators in a magnetic field in a series of 
papers |U E] . Dubrovin and Novikov |H] studied the spectrum of the Pauli operator in 
a two-dimensional lattice with a periodic magnetic field and intensively analyzed the gap 
structure above the ground state. Florek [§J [TU] constructed tensor product representations 
of the MTG to analyze a three-particle system in a lattice in a magnetic field. Kuwabara 
[TT1 IT2] has been studying quantum-classical correspondence from the viewpoint of spectral 
geometry. For example, he ^T] proved that if the whole set of level spacings of the quantum 
spectrum is not dense in R, every trajectory of the corresponding classical particle is a closed 
orbit. Arai [13] found a quantum plane and quantum group structure in the quantum system 
in a singular magnetic field. Thus we can see that quantum mechanics in a magnetic field 
has been an active research area. 

Our study on quantum mechanics in magnetic fields originates from studies of extra- 
dimension models of the space-time. In extra-dimension models the space-time is assumed 
to be a base space of a fiber bundle with a compact fiber or a noncompact fiber. The history 
of extra-dimension models is rather old, but an interest in these models is recently renewed 
as Arkani-Hamed, Dimopoulos, and Dvali ^3] pointed out that the extra-dimension model 
may solve the hierarchy problem of high energy physics. Inspired with extra-dimension 
models we [13] built a model which has a circle S 1 as a fiber over an any- dimensional space- 
times R D ~ 1 . Then we found that a twisted boundary condition in the S^-direction causes 
spontaneous breaking of the translational symmetry. Based on this observation, we [16 
proposed a new mechanism of supersymmetry breaking. Next we |17| ITS] built a model 
which has a two-dimensional sphere S 2 as a fiber over the four-dimensional space-times 
-R 4 . We solved dynamics in the sphere in a magnetic monopole background and then found 
that the monopole induces spontaneous breaking of the rotational symmetry and the CP 
symmetry. We also built a model which has an n-dimensional torus as a fiber and tried to 
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analyze dynamics in the torus in a background magnetic field. However, its analysis was not 
a straightforward task. Then we studied the symmetry structure of quantum mechanics in 
the torus in the magnetic field. We ^J] constructed the MTG in the n-torus and classified 
irreducible representations of the MTG. 

Armed with these tools we are now ready to solve the spectral problem in the n-torus 
T n = R n /A. We decide to solve the problem exhaustively; in our treatment the dimensions 
of the torus is taken to be arbitrary, lengths and angles of edges of the unit cell of A are 
arbitrary, and an arbitrary constant magnetic field is applied to the torus. Thus we aim to 
solve the widest class of quantum mechanics in the n-torus in uniform magnetic fields. 

For busy readers here we give guides for quick access to main results. In Sec. 2 we 
provide a geometric setting to define the problem. The problem to be solved is the eigen- 
value problem of the magnetic Laplacian (|2.13|) with the twist condition (|2.5|) . In Sec. 3 
we find a family of operators that commute with the covariant derivative. Actually they 
are composition of ordinary displacements and gauge transformations as shown at ()3.2|) . 
These displacement vectors form a restricted family of vectors as shown at ()3.8|) . These 
displacement operators generate the magnetic translation group (MTG), which is noncom- 
mutative as shown at (|3.12|) . Along (|3.16|) - (|3.28|) we construct irreducible representations of 
the MTG. In Sec. 4 we introduce a coordinate system, which will be revealed to be useful 
later. In Sec. 5 we define the magnetic algebra by adding differential operators ()5.2|) - ()5.4|) 
and multiplicative operators (|5.9|) to the MTG. Then we construct and classify irreducible 
representations of the magnetic algebra. Sec. 6 is devoted to calculation of simultaneous 
eigenfunctions (j6.8|) of a maximal commutative subalgebra of the magnetic algebra. Then 
we obtain a complete orthonormal set of functions over the magnetic torus, which provide 
an extension of Fourier analysis for the magnetic torus. This is one of main products of 
this paper. In Sec. 7 by applying this method we solve the original problem, the eigenvalue 
problem of the magnetic Laplacian. There we obtain a whole set of eigenfunctions (J7.11j) 
and eigenvalues (j7.12|) . These are main results of this paper. 

2 Gauge field in the torus 

Let t = (t 1 , . . . , t n ) denote a coordinate of an n-dimensional torus T n = R n / Z n . Namely, a 
point (t 1 , . . . , P + 1, . . . , t n ) is identified with (t 1 , . . . , P, . . . , t n ) in T n . A uniform magnetic 
field is generated by the gauge field 

n -r n n 

A = J2 A kdt k = -J2 <t>i* tJdtk + E <*k (2- 1 ) 

fc=l ^ j,k=l k=l 

Here {<pjk = ~<Pkj} and {a?} are real constants. Then the magnetic field is 

1 n 

B = dA = - 4>jk dP A dt k . (2.2) 

2 j,k=i 
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Therefore, the number <ftjk represents magnetic flux which penetrates the (P , t fc )-face of the 
torus. We call the array of numbers (<f>jk) a magnetic flux matrix. 

Let us introduce a complex scalar field / in the torus. The scalar field couples to the 
gauge field via the covariant derivative 

Df = df- 2mAf. (2.3) 

We put the coefficient 2m in front of A for later convenience. Topology of the torus imposes 
a boundary condition on the scalar field. The gauge field itself is not a periodic function on 
R n but it changes its form as 

A{t\...J + l ) ...,t n ) = A{t\...,t\...X) + \fl<t> j kdt k . (2.4) 

1 k=i 

Therefore, if we make the gauge transformation 

f{t\ ..../'• 1 /") e"SL. f^f{t\ ..../' /"). (2.5) 

the covariant derivative ()2.3j) remains covariant as 

Df(t\ . . . , t j + 1, . . . , t n ) = e^SLi +»*Df{t x , ...,t j ,...,t n ). (2.6) 

We call the condition ()2.5|) a twisted periodic condition. There are two ways to bring a 
point (t 1 , . . . , t j + 1, . . . , t k + 1, . . . , t n ) to (t 1 , . . . , t j , . . . , t k , . . . , t n ). The first way is 

f(t\...,P + l,...,t k + l,...,t n ) 
= e *i{4>jk+YZi ^ tl }f(t\ . . . , t j , . . . , t k + 1, . . . , t n ) 

= e «{^+Er=i^i* i +Er= 1 ^* ! }/(t i , ...,t j ,...,t k ,.. . ,r). (2.7) 

The other way is 

f(t\...,t j + l,...,t k + l,...,t n ) 
= e *i{<l>ki+IZ:i ^f{t\ ...,t j + l,...,t k ,...,t n ) 

= ^{^j+Eti^+Eti^f^t 1 , ...,t j ,...,t k ,.. . ,r). (2.8) 

To make these two expressions coincide we need to have 
namely 

e m((f> jk -(f> kj ) _ e 2ixi4> jk _ i (2.9) 

Therefore, compatibility of the periodic conditions ()2.7|) and ()2.8j) demands that <pjk is an 
integer. Hence, the magnetic flux through each face of the torus is quantized. We call the 
torus where the magnetic field has been introduced a magnetic torus. 
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Since two displacements P > P : + 1 and t k > t k + 1 are commutative, we can write the 
twisted periodic condition ()2.5|) in a more general form 

f{t + m) = e m ^l^^ kmHk f(t) (2.10) 

with an arbitrary m = (m 1 ,...,m n ) G Z n . An inner product of two twisted periodic 
functions f(t) and g(t) is defined by 

(f\g}= fdt 1 -- - I' dt n r(t)g(t). (2.11) 
Jo Jo 

Equipped with this inner product the space of twisted periodic functions becomes a Hilbert 
space. 

To define the Laplacian we need to introduce a metric into the torus. Let A be an n- 
dimensional lattice in the Euclidean space R n . We equip the torus T n with a Riemannian 
structure by identifying T n with the quotient space R n /A. Let {ui, . . . , u n } be a set of 
vectors that generates the lattice A. Their inner products is denoted by 

g jk = (uj,u k ) (2.12) 

and its inverse is denoted by g J . Then the magnetic Laplacian is defined as 

A/ = t/" - - 2mA k "j f. (2.13) 

It is also referred as the Bochner Laplacian in literature. The purpose of this paper is to 
solve the eigenvalue problem of the magnetic Laplacian accompanied by the twisted periodic 
condition (|2.5j) . 

3 Magnetic translation group 

Our goal is to find a complete set of eigenvalues and eigenf unctions of the magnetic Laplacian 
(I2.13j) as announced above. A royal road to solving an eigenvalue problem is to detect sym- 
metry. In this section we determine a group of operators that commute with the Laplacian 
and construct irreducible representations of the group. 

The vector space R n acts on the torus as isometries. However, the gauge field restricts 
the admissible class of vectors as seen below. An arbitrary vector v G R n displaces the 
gauge field ()2.1j) as 

A{t) ^ A(t -v)= A(t) - d \\fl <f>jkV j A . (3.1) 

If we perform a gauge transformation of the scalar field simultaneously with the displacement 
f(t) » f'(t) = (U(v)f)(t) = <">^ < ^ k f( t - v), (3.2) 
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then the covariant derivative is changed covariantly 

Df(t) i — > Df{t) = e m ^l^^ kvH \Df)(t - v). (3.3) 

In other words, the transformation U (v) commutes with the covariant derivative as 

(DU(v)f)(t) = (U(v)Df)(t). (3.4) 

Hence it commutes with the magnetic Laplacian, which is defined in terms of the covariant 
derivative. The operator U(v) is unitary with respect to the inner product (|2.11j) . 

The displaced function ()3.2|) also must satisfy the twisted periodic condition. If the 
original function / satisfies the condition ([2.10)1 . the displaced function changes its form as 

f(t + m)=e ^j,fc=i Vjt e ^j,*^^* f (t) (3.5) 

for m G Z n . Thus the displaced function satisfies the condition (|2.1U|) if and only if 

E <f>jkv j m k (3.6) 
j,k=i 

is an integer for an arbitrary m G Z n . In other words, 

n 

J2<f> jk v j (k = l,...,n) (3.7) 

3=1 

must be an integer. We call such a restricted vector v a magnetic shift. The set of magnetic 
shifts forms an Abelian group 

n n = {v e R n | <pv e Z n }. (3.8) 

There is a sequence of Abelian subgroups Z n C Q n C R n . In particular, an integer vector 
m e Z n induces a displacement 

(U(m)f)(t) = e m ^l^ kmHk fit - m). 

However, owing to the twisted periodic condition ()2.10|) . this is reduced to the identity 
transformation 

(U(m)f)(t) = f(t). (3.9) 
Thus we conclude that the group G of effective transformations is generated by 

{U(y) | v e Q n /Z n }. (3.10) 

We call the group G the magnetic translation group (MTG). A product of transformations 
is 

U(v)U(w) = e m ^l^^ kvJwk U(v + w). (3.11) 
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Their commutator is 



U(v)U(w)U(-v)U(-w) = e 2m ^l^ kv3w \ 



(3.12) 



We can say that the MTG is a central extension of the Abelian group Q n / ' Z n by U(l). 

We can express the MTG in a standard form. Let L(n, Z) denote a group of the n- 
dimensional matrices {5*} of integers such that detS = ±1. The matrix S G L(n, Z) acts 
on t G R n by t i — > St and this action induces an automorphism of the torus T n = R n j Z n . 
It also induces a transformation of the magnetic flux matrix as 



E 

i,p=i 



(3.13) 



The Frobenius lemma [20] tells that for any integral antisymmetric matrix there exists a 
transformation to bring it into a standard form 



/ q x 
-qi 



q 2 
-q 2 





-Qn 



qm 





V 



o 



(3.14) 



where {q^} are positive integers and they constitutes a sequence gil^l • • • \q m , which implies 
qi divides qi+\. For example, we may have a sequence 3|6|12|48. Of course, 2m < n. The 
vector subspace of the zero eigenvalue of the matrix <f> has dimensions n — 2m and it is called 
null directions. In the following we suppose that the flux matrix is in the standard form 

dm. 

Now we can write the magnetic shifts (|3.8|) in a more explicit form. Let {ei, . . . , e n } be 
the standard basis of R n in the (t 1 , . . . , t n )-coordinate. Then any magnetic shift is uniquely 
expressed as 



v 



E 



S2j-1 



S2 \ n ~^ m 

e 2j-l H ~ e 2j J + E @2m+k e 2m+k 

) k=i 



(3.15) 



i=i v <U 

with integers {si, . . . , s 2m } and real numbers {9 2m +i, ■ ■ ■ , 9 n }. Namely, the magnetic shifts 
are generated by {(l/^)e 2j -_i, (l/qj)e 2 j \ j = 1, . • • , m} with integral coefficients and {e 2m +fe | k 
l,...,n — 2m} with real coefficients. Hence, if the flux matrix <fi has null directions 
n — 2m > 0, the MTG has a continuous component. Otherwise, the MTG is a completely 
discrete group. 
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Here we summarize our discussion; the MTG is generated by the unitary operators 

U s = U f^-i), V i = U (^) (j = h-..,rn) (3.16) 

and 

W k (6) = U(9e 2m+k ) (k = l,...,n-2m). (3.17) 
According to ()3.9|) . 1)3.11)1 . ()3.12j) and 1)3.14)1 . these generators satisfy the following relations 

(U^ = (V^ = 1, (3.18) 
UjVjU^Vf 1 = e 2ni/q \ (3.19) 
W k (l) = 1, (3.20) 
W k (d)W k (6') = W k (6 + 6') (3.21) 

and other trivial commutators. 

To solve the eigenvalue problem of the magnetic Laplacian we need to prepare the whole 
set of irreducible representations of the MTG. Let {|r 1; . . . , r m ; dx, . . . , d n _ 2m )} be elements 
of a representation space that are labeled by 

r.eZ/Z^ (j = l,...,m), (3.22) 
d k eZ (k — I, . . . ,n — 2m). (3.23) 

Then the generators (j3.16j) and (j3.17j) are represented by 

Uj\n, . . . ,rj, . . . ,r m ; dx, ■ ■ .,d n - 2 m) = \rx, ■ ■ ■ , rj + 1, . . . , r m ; dx, ■ ■ .,d n - 2 m), (3.24) 
Vj|rx, . . . , Tj, . . . , r m ] di, . . . , d n —2m) = e 3 ^ l^ij ■ ■ ■ ■> fj, ■ ■ ■ , f m ] di, . . . , d n —2m)i (3.25) 
W k (9)\rx, . . . ,r j} . . . ,r m ;dx, ■ ■ ■ , d n - 2m ) = e~ 2mdk6 \rx, . . . , r j} . . . , r m ; dx, ■ ■ ■ , d n -2m)- (3-26) 

Thus 

<?i q m 

H d = • • • ® C|n, . . . , r w ; d l9 . . . , d n _ 2m ) (3.27) 

ri=0 r m =0 

provides an irreducible representation space of the MTG. Its dimension is 

d\mH d = qx X g 2 x ■ ■ • x q m . (3.28) 

The labels rj and d k in ()3.22)) and ()3.23)) will become good quantum numbers for the 
Laplacian (j2.13j) . The dimension ()3.28)) will give the degree of degeneracy of each eigenvalue. 

4 Diagonalization of the magnetic field strength 

To solve the eigenvalue problem of the magnetic Laplacian we need to take the metric ()2.12)) 
into account. Remember that the basis {ux, ■ ■ ■ , u n } generates the lattice A in the Euclidean 
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space R n and that the torus is isometric to R n / A. Let x = (x 1 , . . . , x n ) be an orthonormal 
coordinate of R n . It is related to the normalized coordinate t = (t 1 , . . . , t n ) via 



x = t l u x H h t n u n = Ut. 

In these coordinates the magnetic field (|2.2|) is expressed as 



(4.1) 



B = 2 E B Jk dx j Adx k = - £ B 3k A u k p dt l Adt p = -J2 K dt l A dt p . (4.2) 

j,k=l j,k,l,p=l l,p=l 

The number Bjk is areal density of magnetic flux which penetrates the {x\ x fc )-plane. If we 
perform a coordinate transformation 

x = Ry (4.3) 

by an orthogonal transformation R G 0(n,R), the components of the field strength is 
transformed as 



1 71 1 
B = - Bjk dx j A dx k = - 

j,k=l j,k,l,p=l 



£ B 3k R\R\dy l Ady?. 



(4.4) 



By a suitable orthogonal transformation the field strength matrix (Bjk) can be brought into 
a standard form 

{ vi \ 
-v x 

v 2 
-v 2 



v = RBR 



z/ m 
-i/ m 







(4.5) 



where {yj\ are positive real numbers. The (i 1 , . . . , t n )-coordinate system block-diagonalizes 
the magnetic flux in the form of ()3.14|) while (y 1 , . . . ,y n ) block-diagonalizes the magnetic 
field strength in the form of ([4.5)1 . The transformations (|4.1|) and (j4.3|) are combined into 



y = R X = RUt = Lt. 
Then we obtain relations among the matrices 

cj) = HJBV = HJRv RU = %pL. 
The phase factor in (|3.2|) is rewritten as 



(4.6) 



(4.7) 



j,fc=l j,k,l,p=l j,l,P=l 3=1 '=1 

(4.s; 



ii 



The gauge field (|2.1|) is expressed in the ^-coordinate as 



3=1 



E «i(£ _1 )W 



We put 
for later use. 

Actually, we can choose a transformation matrix L that has zeros in this pattern 

/ r1i-\ i-2i-l 

2p-l ^ 2g 1 

L - 



(4.9) 



(4.10) 



^ 2p-l ^ 2g 
r2i r 2j 



J_/ 2m+r 
Lj < 



2p-l ^ 2q ^ 2m+r 

j 2m+k t 2m+k j 2m+k 
\ 2p-l 2g 2m+r 





f * 





\ 


)- 


* 


* 







K * 


* 


* J 



(4.11) 



with i, j, p, q — 1, . . . , m and k,r = 1, . . . , n — 2m. The inverse matrix L~ l also has the same 
pattern of zeros. This distribution of zeros is proved in the appendix A. Then we can rewrite 
the magnetic shift operators (|3.16|) and (|3.17|) with a help of (|4.8|) in the y-coordinate as 

'21-1 ,,,,21 j2l ,, „,2l-l\ „ / „■ „• , 

(4.12) 

(4.13) 
(4.14) 



Vjftf) = e--^^)n-^- 1 /( ^_ L i 2 . (1/g . ))) 

wmftf) = f(y l -L l 2m+k 6) 



for j 



m and k 



v % l, 2j _ 1 



n — 2m. From (|4.11j) and (J4.7|) we get a formula 

m 

2^ V l L 2j-l L 2p\ L ) 2i 
,p=l 
m 



-l\2p 
2i 



l,p=l 
m 

p=l 



-l\2p 
) 2i 



(4.15) 



which will be repeatedly used later. 



5 Magnetic algebra 

In this section we introduce new operators which act on twisted periodic functions. These 
new operators and the operators in the MTG generate an algebra, which we call a magnetic 
algebra. We construct and classify its irreducible representations. In the next section we will 
prove that the space of twisted periodic functions is actually an irreducible representation 
of the magnetic algebra. In this sense, the magnetic algebra completely characterizes the 
quantum mechanics in the magnetic torus. 
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Now we introduce a family of Hermite operators. Expanding the covariant derivative 
(12. 3 j) in terms of the ^-coordinate 



Df = iY,Pifdy l , (5.1) 
i=i 



we define differential operators 



P23-1 = -i[^- T + 7TW j y 2j -27ri^ 1 ) (5.2) 
P 2j = -i (-^j - mujy 2 ^ 1 - 2irif3 2j \ (j = 1, . . . , m) (5.3) 



D 
dy 2 



^2m+k = -t I » 2m+fc - 27rz/3 2m +fc (A; = 1, . . . , n - 2m). (5.4) 



These are Hermitian with respect to the inner product (|2.11j) . Since (y , . . . , y n ) is an 
orthonormal coordinate, the Laplacian (j2.13j) becomes 



- Af = E(^) 2 /- (5.5) 
i=i 

Nontrivial commutators among P's are 

[P 2j ^,P 2j ] = 2mu J (j = l,...,m). (5.6) 

The other commutators vanish. We call the operators {P 2 j-\, P 2 j} transverse momenta 
while we call the operators {P 2m +k} longitudinal momenta. Since the covariant derivative 
commutes with the magnetic shifts as seen at (|3.4|) . the momentum operators {Pi} commute 
with the shift operators {Uj, Vj, Wk}, which are defined at (|3.16|) and (|3.17|) . 

Next we introduce another family of unitary operators. For this purpose we need an 
observation; in the t-coordinate that expresses the magnetic flux matrix in the standard 
form (j3.14|) . t 2m+h (k = 1, . . . ,n — 2m) are genuine cyclic coordinates. That is to say, the 
twisted periodic function (|2.5j) is periodic with respect to these coordinates as 

f{t\ . . .,t 2m+k + 1, . . . ,t n ) = f(t\ . . .,t 2m+ \ ...,t n ) (k = 1, . . . ,n - 2m) (5.7) 

and the magnetic shift (J3.2j) is reduced to an ordinary continuous shift 

(W k (9)f)(t\ . . . , t 2m+k , ...,t n ) = f(t\ . . . , t 2m+k -6,...,t n ). (5.8) 

Then we define an operator T k for each k — 1, . . . , n — 2m which acts on / by multiplication 

(T k f)(t) = e 2mt ' 2m+k f(t) = e^t^-^y'fit). (5.9) 
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Here we used the inverse L 1 of the coordinate transformation ([4.6)1 . The operators {T k } 
are unitary operators with respect to the inner product (|2.11|) . They satisfy 

W k (6)T k = e~ 2md T k W k {Q), (5.10) 
[P i ,T k } = 2n(L- 1 ) 2 t m+k T k (i = l,...,n; k = l,...,n-2m) (5.11) 

and commute with the other generators of the MTG. 

Combining all the operators introduced above we define an algebra A with the generators 
{Pi, Uj, Vj, W k {6),T k \i = 1, . . . , n; j = 1, . . . , m; k = 1, . . . , n - 2m; 6 e R} and with 
the relations lETTS|) - lEr2T|) . (pTKj) . lETTTO) . (fSHH) and other trivial commutators. We call 
the algebra A a magnetic algebra. In the following we will construct all the irreducible 
representations of the algebra A and classify their unitary equivalence classes. 

A subset of generators {P 2 j-i, P2m+i, Vj, Wi{6) \ j = 1, . . . , m; / = 1, . . . , n — 2m; 9 6 R} 
generates a maximal Abelian subalgebra of A. Hence these generators are simultaneously 
diagonalizable. Their simultaneous eigenstate \k,r,d) is labeled by r 3 - G ZjZ q . of ()3.22|) 
and d\ G Z of ([3.23)1 with new labels 

k 2 j-i,k 2m+ i e R. (5.12) 

The generators {P2J-1, P2m+i} act on these states as 

P 2 j-i\k,r,d) = 2nk 2 j-i\k,r,d), (5.13) 
P 2m +i\k,r,d) = 2iik 2m+l \k,r,d), (5.14) 

and {Vj,Wi{6)} act as ([3.25)1 . ()3.26)1 . respectively. The coefficient 2n was put for later 
convenience. Other generators {P 2 j,T l | j = 1, . . . , m; I = 1, . . . , n — 2m} act on the states 
as 

d 

(ij\P 2j \k,r,d) =ivj-=L (tp\k,r,d), (5.15) 

OK 2 j-i 

T%, r, d k ) = e 2-{E7 =1 (^-i+A i fc 2 ,_ 1 )A'fc 2 ,M + E::r z lp dp}\ k . + A i kh r? dk + ( 5 . 16 ) 

with A l ki = (L~ 1 ) 2 J n+l . We will determine the (n — 2m) x (n — 2m) matrix Z lp soon 
later. Here represents an arbitrary state. The rests {Uj \j = 1, . . . ,m} act as ()3.24[) . 
Therefore, in an irreducible representation space the eigenvalues k 2m+ i are linked to di via 

n— 2m 

k 2m +i = d>t(L 1 ) 2 2 ™ 4 fi — P 2m +i 

1=1 
n—2m 

= E ( d i- «2m+i)(^ _1 ) 2 2m+i (i = 1, - - - , 71 - 2m). (5.17) 

Z=l 

Here the real number (3 2m+i coincides with the one that appeared in ()5.4j) . We also used the 
relations P~TU)1 and (|4~TT)l. 
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The matrix Z lp is determined by the condition T l 'T l = T l T v '. The action of T l 'T l gives 
T l 'T l \ki, r, d k ) = e ^{££i(^- 1+A '^- 1 ) A W^+££: 1 2m z ' Pfi W 

e 2^{^ j m =1 (fc 2j _ 1 +A i fe 2j _ 1 +A''fc 2:) _ 1 )A i 'fc 2 ,/^+^;: i 2m Z''r(d p +5j)} 

\ki + A% + A 1 ' h, r, d k + 5 k l + 5 k 1 ') (5.18) 
while the action of T l T 1 ' gives 

T l T l '\ki, r, d k ) = e 2 " 1 ^^^'^''^-'^''^'^^!^^''^} 

e 2«{^; 1 (fc 2j _ 1+ A^ 2rl +A'fe rl )A ! t 23 /, J +E;: 1 2m z'f(dp+i/)} 

+ A''fci + A'fci, r, 4 + 5/ + 5 k l ). (5.19) 
To give T ? 'T' = T l T v the matrix Z"' must satisfy 

m m 

Y,(&%j-iA!'k2j/vj) + Z 1 ' 1 = Y,(A l 'k 2j - 1 A% j /v j ) + Z 11 '. (5.20) 

3=1 3=1 

A general solution of the above equation is 

m 

Z lV = ^2(A l k 2 j-iA 1 ' h 2 j /vj) + S 11 '. (5.21) 

3=1 

Here we leave an arbitrary symmetric matrix S 11 ' = S 1 ' 1 yet undetermined. Actually any 
choice of S 11 ' results in an equivalent representation, and therefore we take S 11 ' = 0. 

A restricted set of vectors {\k, r, d)} that are labeled by the mutually independent pa- 
rameters 

fc 2j _! G R, Tj G Z/Z qj , d t eZ (j = 1, . . . , m; I = 1, . . . , n - 2m) (5.22) 

spans a Hilbert space TC a for each fixed value of (a 2m +i, «2m+2, • • • , ot n )- Thus we conclude 
that a unitary equivalence class of irreducible representations of the algebra A has one-to-one 
correspondence with the parameter (a 2m+ i, a2m+2, ■ ■ ■ ,a n 

6 Fourier analysis for the magnetic torus 

Now let us turn to the space of twisted periodic functions. It is a representation space of 
the magnetic algebra. We will calculate the whole family of eigenfunctions of the maxi- 
mal Abelian subalgebra of the magnetic algebra. These eigenfunctions Xk^dit 1 , . . . ,t n ) = 
{t\k, r, d) satisfy 

P23-I Xk,r,d = ~i (^ Q y 2j~l + KWjV 23 - 2nif3 2 j-l) Xk,r,d = 2ixk 2j -l Xk,r,d (6.1) 
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( d ■_ \ d 

P23 Xk, r ,d = ~i yoyTj - nivjy 2 ^ 1 - 2iri/3 2j \ Xk,r,d = Wj ^ ~ Xk, r ,d (6.2) 

P2m+l Xk,r,d ~ ~i I Qy2m+l ~ ^ifam+l J Xk,r,d 

(n—2m \ 
£ d^L- 1 ) 2 ^ - (3 2m+l \ X k,r,d (6.3) 

Uj Xk, r ,d = e %it23 Xk,r,d [t 2 ^ 1 ~ = Xk,(r+i),d (6.4) 

Vj Xk, r ,d = e- mt23 -\ Kr4 (t* -^j= e- 2mr > /q >Xk, r ,d (6.5) 

W0) X k,r,d = Xk,r,d(t 2m+l ~0) = e~ 2mdl0 Xk , r , d (6.6) 
T l Xk,r,d = e 2 ^ +l Xk , r4 = e-S" 1 (^- 1+ A fc ,_ l) A fc ,/, J+E ;:r ^ Xk+AkMd+1) (6.7) 

for j = 1, . . . , m; I = 1, . . . , n — 2m, and 9 E R. Here (r + 1) is abbreviation of (r; + Sij) and 
(d + 1) is abbreviation of (d p + <5 Z ). And Aki = (L~ 1 ) 2 ™ 1 ' +1 as given in f)5.16j) . Then Xk,r,d is 
a simultaneous eigenfunction of {P 2j -i, P2m+u Vj, Wi(9)}. In the rest of this section we will 
solve the set of equations (j6.1|) - ()6.7|) to get the solutions 

Xk,r,d(y > • • • > 2/ ) = ce ^= lZjH u ^2,-1^2,/ 3 

oo 

E 

(Tl,lJ2,...,0' m =— OO 



e ^'j 
e 



e 2 «Er 1 =i(«^^-^)( i " i )^^-Er=i i "-i(^.+''.)/*} 



e 



_2tt 



(6.* 



where c is a common normalization constant. The coefficients jj are given later at (|6.13|) . 
The fact that the eigenfunctions are uniquely determined up to the common coefficient c 
implies that the space of twisted periodic function is an irreducible representation space 
of the magnetic algebra. Consequently, the eigenfunctions ()6.8|) constitute a complete or- 
thonormal set of the space of twisted periodic functions over the torus. This is one of main 
results of this paper. Hence an arbitrary twisted periodic function over the torus can be 
expanded as 

f(y\ . . . , y n ) = £ X k , r ,d Xk, r ,d(y\ ■■■^ n ) (6-9) 

k,r,d 

with unique coefficients \ k , r ,d- Therefore, the complete set {Xk,r,d} provides a new basis for 
Fourier analysis in the magnetic torus. 
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In the rest of this section we give detailed lengthy calculations to prove the above state- 
ments. The reader may skip them to the next section, where we calculate solutions of 
the eigenvalue problem of the magnetic Laplacian using the main result (|6.8j) . First, a 
simultaneous solution of ()6.1j) and ()6. 3j) is 



^ ^'^.Xv 2 , y\-.., v 2m ), (6.10) 



where 4>k,r,d(y 2 , 2/ 4 , ■ ■ • , y 2m ) is an arbitrary function to be specified later. 

Next let us turn to the other equation ()6.5|) . Using ()4.13|) we can rewrite ()6.5|) as 

e-rtE^O/^^ _ L * 2 •/ g^) = e- Sh<r ^«Xik,r 1 d(y < ). (6-11) 



As discussed at (j4.11|) we have taken the matrix L such that L 2 j =0. Therefore, when 
(|6.1U|) is substituted, the LHS of (|6.11|) becomes 



>Er= 1 {fe- 1 +ft i -i)!' 2, - i -(V2)^ 2i - i i/ 2 %^^r^(i- i )t + + p 1 !/ 2 "' +! 



2jlH 3 , 
2 

e 



o.' V™" 2 " 1 J f t —l\ 2m +P T 2m + l / „. „. 

e -^E P ^ML ) 2m+l L 2j l^ (j)k ^_ L y q .y (6 .i2) 

Hence, if we put 

n-2m 

7 , = £ d p (L-^L 2 %+ 1 , (6.13) 
p,/=i 

(j6.11|) implies that 

If we introduce another coordinate system (z 1 , z 2 , . . . , z m ) which is related to (y 2 , y A , . . . , ?/ 2m ) 
via 

m 

y 2i = J2L%(zi/ qj ), (6.15) 

then (J6.14j) is rewritten as 

</> k , r , d {z\ . . • , J - 1, . . . , z m ) = e 2m ^~ r ^ </> k , r , d {z\ ...,z\...,z m ). (6.16) 
Moreover, if we put 

^ r4 {z\ ...,z m ) = e 2 «Er =1 (7i-^>V^ ^ krd{z ^ . . . f z m^ (6 17) 

then (|6.16j) implies that 

• • • , ^' - 1, • • • ,^ m ) = Vw(^\ . . . , . . . , z m ). (6.18) 
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Hence ipk,r,d is a periodic function with the period 1 and can be expanded in a Fourier series 

oo 

^ Kr4 {z\ z m ) = J2 <W e 2 ™ Ej=1 (6-19) 

Note that the inverse transformation of 1)6.15)1 is given by 



m 



[z 



/Qj)=J2( L )iiV 2i - ( 6 -20) 
i=i 

Combining the above equations we can write down the eigenfunction (JbMUj) in a more specific 
form as 

g Z_-/j,Z — 1 J x ' 2m+I y 

E Cfc^e^^r^ife^^-^)^" 1 )^^. (6.21) 

(Tl,(T2,...,0- m = — CXD 

Moreover, referring to ()4.11|) and (|4.14jl . we can see that ()6.21|) satisfies ()6.6|) . Thus we have 
seen that Xk,r,d is a simultaneous eigenfunction of {P 2 j-i, P2m+u Vj, Wi(9)} as announced 
above. 

The remaining task is to solve (|fi.2j) . ()6.4j) . and (jfi.7jl . Let us begin with (|fj.2j) . The LHS 
of (EH is 



r,c/ 



<ti,iT2)...,0'ot=— oo \ i=l / 



e 

21 



u fe,r,a,<T ° 

On the other hand, the RHS of is 

in-?— & r d = in eW^-^^-a^w^E^ 



O r l,0'2,---iO"m = — OO 

Therefore, we have an equation 

= 2?r g( ft(T . + r . _ 7l )(L- 1 )%. - ) c fc , Mi<T (6.22) 

0*23-1 \i=l ' ) 

and get its solution 

C k ,rA° = C ^Aa e ~ 2 ^E™ ^fe^-^" 1 ^ k^-lM-E^l ftiM'rf ( 6 . 23 ) 
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Thus (I6.21J1 becomes 

Xk rd(y\ y 2 ,---, y n ) = e^^Ai^^+^-^^+k^i^/^^m^^y 23 } 



n ~ 2m rJ .( T -l\2m+j 2m+i 



oc 



E c 0<r , d , a e m ^=^ aj+rj -^ ){L )il(y - k ^l^\ (6.24) 



^u,r,a,<T * 

CTl,CT2,...,(T m = — OO 

Next we turn to (fOjl . With the aid of (l4~T2l) and (ETToT) the LHS of JS3|) becomes 



= ^e^i^-i^v^-^xl- 1 ^ 

oo 

V r„ ^ p -2^Eri=ife CTl+ri -^)( L " 1 ) 2 2i L2 2 J -i/^ 

<Jl,(T2,...,cr m =— oo 

To make this coincide with the RHS of ()6.4j) we have a recursive equation 



e 



it is symmetric as 



(6.25) 



e-^E^xto^-T*)^- 1 ^-!/* COirjdi(j = Co , (r+1)idiCT . (6.26) 
Here (r + 1) means (r* + If we define an m x m matrix 

m m 

^ = E(^" 1 ) W-i/* = E ^ ar-iiVi/(»»). ( 6 - 27 ) 



- = E( L «"-i^^Vi - L\-& L %±i)Kw) = <k*-w-i/(<kQs) = (6-28) 
i=i 

by virtue of (|4.7jl . Then the solution of ()6.26|) is 



c 0,r,d,cr — CQ,0,d,CT e 

e~ 

e -27riX;™, i= ife^ + (V2)n-7 I )(i- 1 ) 2 2 I ; L 2 ^_ 1 »- J /g J 

Therefore, ()6.24|1 becomes 

Xkrdiy 1 ^ 2 , ■■■,y n ) = e^E^JC^-i+^-O/^'+^-i^M-Ci/a)^-^'- 1 ^'} 

e -2«Er^i fe-ii^r,/^ e «E™, i= i ^^K-i^-ifri/fli)^/*) 
oo 

E C 0,0,d,CT 

>r — 1 \2i / 9/ . / \~^m T 97 / \ 

(6.30) 



<Tl,<T2,...,Cr m =— oo 
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Since (Uj) qj = 1 by ([3. 18)1 . the substitution rj \— > rj + qj must leave Xk,r,d invariant. This 
substitution gives 

Xk (r+q) d{y) = e 2 "S-i {(fc2l - 1+/32l ' l)y2 " 1+fc2l - l/32lM_{1/2) ^ y2l " ls/2l} 

e - 2 ™^Z=i 02i-iL 2 ;-_\ ri / qi e «E™ p=1 -nt\ £ V( p «/«)( r »/*) 

OO 

^ C 0,0,d,(o--1) 



(Tl,(T2,...,0" m = — OO 



e -2 7 riX;™ i=1 fe^+n-7 l )(l/9 l )^i2 l " 1 1 i 2 ^_ 1 

Here (cr — 1) is an abbreviation of (crj — %). We used f!4.15|) . Then we have another recursive 
equation 

e -27ri^™ =1 (g lCTi +r 1 -7 i )(l/g0^i 2 2r_ 1 i i2 2j-i (6 32) 



Remember that qiqjY^ = E/=i ^^i-i-^y-i = QjQi^ji is symmetric. The solution of (|6.32j) 



is 



Co,0,d,<7 = Conine Z-,,,j=i Z-,,=i »V ^ 2™+; 2,-1^ Z^, J= i^i 2,-1. 



e -27riEr^ 3 = l((V 2 )9 iCTi+r »-^)( 1 / g »)^ i2 2 l ~A L2 2i-l <T J (g 33) 

Substituting it into (j6.30|) and using ()4.15|) . we get 

CTl,(T2,---! CT m = — OO 

g 7 ™ Ejj,j=i v l L %i-\ L %-\ (<li<7i+ r i)(qjVj+rj)/(qiqj) 

e^E^ito^^-^J^'^M^-E^x^K-i^-H-iJ/ft} 



g t—/i,l=l — \ ' Zm+I i» -i— 'J = 

g 2 ™ E Jli fc 2 j - 1 {j/ 2j ' " 1 +/?2j /Vj ~YX= 1 ^M-l (« CT ' _ Ti )/« } (6 34) 



20 



Finally, we are going to solve (|6.7j) . Its LHS becomes 

oo 

E c 0,0,d,0 

<Tl,<T2,...,Cr m =— oo 

EH »=1 U P L 2i-l i2 27-l(* 'i+ T 'i)(9j a 'j+ r j)/fe9j) 



2 

e 



e ^p> ; 

e^E;; 2 ^^- 1 )^^^ 

■\n — 2m v — \m 

e 



2^Er=i( fc 2i-i+(i- i ) 2 2 7^)te 2j - i +/?2 J /^-Er=i L %"- i i(9p^+^-7p)/9p} 

-27ri^^ 1 (L- 1 ) 2 ™+{{/32 J /^-^™=i L2 2VA(9p CT P+ r p->)/9p} (g 35) 



We put A z fcj = (L i) 2 ™" 4- ' as before. The change d p ^ d p + 5 p l causes a change of 7j as 

n— 2m 



^ = E (l- 1 ) 2 ^ 2 ^ 

n m m 

= E^ 1 )^^ - E(£ -1 )^V - E( L_1 ) 2 2T + ^%- 

i=l i=l i=l 

m 

= - - J2(L-^ +l L% (6.36) 

i=l 

via (|6.13|) with (|4.11j) . Moreover, we can see that 

m m 

E^%-iAV% = - e ^(^nr^v^ 

p=l p,«=l 

m 

= - e (^ 1 )i(^- 1 ) 2 2 r^y^ 
p,i=i 

= -E(^T 2 ? + '4M 
i=i 

= -A l k 2j /vj. (6.37) 



Therefore (I6.35j) becomes 



e 



-^E^i^-^io^M+E^i^-Wsi.) 



21 



oc 



^ Co,o,d,0 



cri,<T2,...,cr m = — oo 

J27j P =l ^ L «- 1 ^fj - 1 °"i + r i ) (<?j Cj +rj)/(,cn qj ) 



e 



2^E-;=r^+ 5 D(i- l ) 2 2 ™ + ; P to 2m+p -Er= 1 i 2 2 7-rfe^+^)/^} 



In the course of calculation we used the fact Y%=i{L 1 ) 2 ™ + ^ P 2i-i = 0- With the aid of 
(14.15)) . ()6.36|) and the definition (|5.21j) we can deduce that 

m n—2m 

£ {L-^[L%\ lv lq v = - £ Z li di. (6.39) 

j,p=l i=l 

Thus we reach 

6 Xk,r,d — e J Co,o,d,0 C 0,0,(d+1),0 

e^^-i+A^-OA^^^^ (64Q) 

To satisfy ()6.7|) we meet another recursive equation 

co,o,(d+i),o = e ^= lV ;2j - lP 2j/ J c 0M0 (6.41) 

and we get the solution 

c„,o,,,o = (6.42) 
Substituting it into (|6.34j) we reach the final result 



X k ,rM = 00,0,0,0 e- 2 "^^ * (L 



e - 7 ™E J "li^ 2j V J 



E 



(Tl,(T2v! CT m = — OO 



e 

'77/ > 

e 



(6.43) 

This is the result (16. 8|) announced previously. The eigenfunction Xfc,r,d is determined up to a 
unique normalization constant Co,o,o,o- Thus we conclude that the space of twisted periodic 
functions over the torus is an irreducible representation space of the magnetic algebra A. 
Finally, we have proved that the set of functions {Xk,r,d} is a complete orthonormal set in 
the space of twisted periodic functions over the torus as announced at (|6.9|) . 
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7 Eigenfunctions of the magnetic Laplacian 




In this section we will write down explicitly solutions of the eigenvalue problem of the 
magnetic Laplacian (J2.13j) or (J5.5j) . which is expressed in the ^-coordinate as 

a/ = U I + ™>y 2j - 2 ™/Vi) / + - ti^- 1 - 2™/%) / 

2m+fc - 2 ™/w) /• (7-1) 

The function / must satisfy the twisted periodic condition (|2.5jl . We will obtain solutions 
using the Fourier analysis that is developed in the previous section. 

Here we would like to describe the outline of our method. As mentioned in the previous 
section, the Laplacian commutes with the magnetic shift operators, {Uj, Vj, Wk\- Hence 
the labels (r, d) = (r±, r 2 , . . . , r m , di, d 2 , ■ ■ ■ , d n _ 2m ), which are defined in (|3.22j) - (|3.26|) . are 
'good' quantum numbers. Moreover, the Laplacian commutes with the longitudinal momen- 
tum operators {P 2m +k}- Hence the corresponding momentum eigenvalues {k 2m +k} are also 
good quantum numbers and are related to the labels d\ via (j5.17|) . On the other hand, the 
Laplacian does not commute with the transverse momentum operators {P 2 j~i, P 2 j}- Hence 
the transverse momentum eigenvalues {&2j-i} do not remain good quantum numbers. The 
Laplacian admits a new set of good quantum numbers (ni,n 2 , ■ ■ ■ ,n m ), which will be in- 
troduced soon later. It will be revealed that eigenfunctions of the Laplacian are actually 
matrix elements of a unitary transformation, 

ip n ,r,d{ k ii • • • > k 2m -i) = (k,r,d\n,r,d), (7.2) 

which relates the quantum numbers n's to fc's. In the /c-space it is rather easy to get 
eigenfunctions by the standard method of a harmonic oscillator. On the other hand, the set 
of eigenfunctions of the momenta and magnetic shifts, 

Xk,r4iv l i y 2 ,---, y n ) = (y\k, r, d), (7.3) 

plays a role a unitary transformation which bridges between the momentum space and the 
real space like the usual Fourier transformation. Hence the Laplacian eigenfunctions are 
transformed into the ^/-coordinate representations by 

rod 

^n,r4(y l iV 2 i ■ ■ ■ iV n ) = (v\n,r,d) = / dkxdk^- ■ ■ dk 2m -i(y\k,r,d) (k,r,d\n,r,d). (7.4) 



This will give the desired result. 

Now let us carry out the program outlined above. We define creation and annihilation 
operators associated with the transverse momenta as 

fljt = ^i^ 2 ^ 1 ~ iP2j ^ aj = ti^^ 1 + 0' = • • • > m )- ( 7 - 5 ) 



23 



It is easily verified that [aj,a£] = 6jk- Then the Laplacian ()5.5|) becomes 



n— 2m 



(7.6) 



The eigenstate for the lowest eigenvalue satisfies 

1 



fc=l 



= (fc|a i |/2) = - 7 =(A;|(P^_i+iP2 J -)l fi > 







^- 1 + Vi dh^, 



(k\n). 



Here we used (|5.1Hj) and (|5.15j) . The solution is 

States for higher eigenvalues are generated by creation operators as 
(k\n) 



(7.7) 



(7.1 



. (k\(a\r-.-(alr-\f2) 
Vni! • • -n m \ 

■y m 



n 



0*1 ! • • - Wm! 7 =1 Lv 75 



1 



2vrA;2j-i - z^- 



<9 



e TE7 =1 (^-o 2 / 



(A;|J7) 
9 



- 27 rE7 =1 (^-i) 2 A 



(7.9) 



for ni, ra 2 , . . . , n TO = 0, 1, 2, . . .. We are suppressing other labels (r, d). In the appendix B we 
prove that 



dz 



= 07 (iV2^Y e-™ z2 H n (zV2w). (7.10) 

In the second line H n (£) is the n-th Hermite polynomial. Substituting ()6.8|) and ()7.9|) into 
f!7.4|) and applying ()7.10|) we obtain 

/•oo 

{y\n,r,d} = / dk x dk 3 ■ ■ ■ d 2m -i (y\k, r, d)(k, r, d\n, r, d) 



E 



ffl,(T2,... 1 ffm = -<» 

.. r 2i— 1 



e vi J2Z,i^i^ L2 L-i L %-i(^ a i +r ^(^ a J +r j)/(^<ij) 
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n 



/i2. _ e-^ife^-'+AiM-E^a i ! £l 1 i(»o-i+ri-70/«} a 

(m 
i=i 



•(7.11) 



This is the main result of this paper. Using (|5.17j) . we can calculate eigenvalues of the 
Laplacian ()7.6|) as 



1 



m ( 1 

3=1 V Z 



i n— 2m 

z fc=i 



-l\2m+Z 
J 2m+fc 



tpn. 



r.d- 



n—2m 

Y ^(di - a 2m +i){L~ 
i=i 

(712) 

Eigenvalues depend on quantum numbers n\,n 2 , . . . , n m = 0, 1, 2, . . . and d\, d 2 , ■ ■ ■ , 4-2m 
= 0, ±1, ±2, . . . but not on r = (r 1; r 2 , . . . , r m ) G x Z ?2 x ■ • • x Z 9m . Thus each eigenvalue 
is degenerated by qiq 2 ■ • • q m folds as predicted in (J3.28J) . If a ratio Vijvj{i ^ j) is rational, 
degeneracy happens more. On the other hand, if a 2m+ i = 0, the eigenvalue for — di coincides 
with the one for d\. If a 2m+ i = 1/2, the eigenvalue for (— d\ + 1) coincides with the one for 
di. Moreover, for specific values of {{L~ 1 ) 2 2 n ^ k } we may meet more multi-fold degeneracy. 

Let us discuss physical meanings of the eigenvalue ([7. 12)1 . It is energy of an electrically- 
charged particle moving in the magnetic field in the torus. In 1)7.12)1 we put the coefficient 
—1/2 in front of A to adjust the equation to the conventional Schrodinger equation. In 
the context of classical mechanics, the particle exhibits a cyclic motion with the frequency 
Uj in the (y 2j_1 , y 2j )-plane for each j = 1,2, ... ,m. And it exhibits a uniform straight 
motion along the y 2m+k -ax\s. The whole motion is a superposition of those cyclic and 
straight motions. When we turn to quantum mechanics, energy of the cyclic motion is 
quantized and results in the so-called Landau level 2irvj(rij + 1/2). On the other hand, 
the longitudinal momentum P 2m +k associated with the straight motion is quantized to be 
2-Kk 2m+k = 2nE?=i m (di-a 2m+l )(L- 1 ) 2 Z + l k = 2tt ^~ 2m di(L~ 1 ) ! ^ k - f3 2m+k with integers 
(di, d 2 , . . . , <i n _ 2m ) as explained in 1)5.17)1 . Along the course of the straight motion the 
particle flies around the torus and picks up the so-called Aharonov-Bohm effect. Then the 
momentum is shifted by the Aharonov-Bohm parameters {(5 2m+ \, . . . ,(3 n ). Accordingly, the 
kinetic energy of the straight motion is also quantized. The total energy is then given as 

(JOT- 

Moreover, let us examine meanings of other Aharonov-Bohm parameters (/3i, . . . ,f3 2m ). 
These do not affect the energy ()7.12)) and hence they have a geometric significance rather 
than a physical significance. To understand their meaning we rewrite the eigenfunction 
ijOSI) as 



Xk,r,d(y) = ce 



- 27 ™Ei,;=i Ay-i£ w-iTi/?i p -2™ E J= i Ez=i 
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oo 



E 



<7l,<72,...,(7 m =— oo 

6^^^3,1=1 ^ i2 i~- 1 i- L2 2 J -l(^°"i+ r i)( , 3j "j+ r 'j)/(9i9j) 



From ()6.36|) we have 

m n-2m m 



Then (j7.13|) is rewritten as 



Xk,r,d\y) = Ce ^ i=1 ^3=1 ,1A y 2 J M 3 y )2 3 -lWS/ 3 



— -m 

e 



'E,"Li^^ 2j V 



(7.13) 



E /Vi^-W?* = - E E (7-i4) 



oo 

E 

(71,(72, ...,<7m = — 00 

g^E^.i^i^^^rA^i-^^^+^oc^^+^^/c^^) 
6 2 -Em=7 ^-^V^-EI^ ^(«"*-H-i)/«} 

Thus we can see that the parameters (Pi, ... , ft m ) induce a displacement 

(yV-\ y») -> (y^ 1 + ^-M, Z/ 2J - /Vi/"i) (7-16) 

of the profile |xfc,r,d(2/)|- This is the geometric significance of the transverse Aharonov-Bohm 
parameters. 



(7.15) 



8 Conclusion 

Here we summarize our discussions. As well-known, a U(l) gauge field replaces the partial 
derivative by the covariant derivative and generates a magnetic field. As a natural extension 
of the eigenvalue problem of the ordinary Laplacian in the n-torus, we formulated the 
eigenvalue problem of the magnetic Laplacian. The ordinary Laplacian admits continuous 
Abelian symmetry and therefore the usual Fourier analysis is applicable. However, the 
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magnetic Laplacian does not admit continuous Abelian symmetry and therefore the usual 
Fourier analysis is not applicable to it. Hence, we developed an alternative method, which 
became an extension of the usual Fourier analysis. We identified symmetry structure of the 
magnetic Laplacian and defined the magnetic translation group (MTG), which is discrete 
and non-Abelian in general. Moreover, we defined the magnetic algebra by extending the 
MTG. We proved that the space of functions on which the magnetic Laplacian acts is an 
irreducible representation space of the magnetic algebra. By diagonalizing the maximal 
Abelian subalgebra of the magnetic algebra we obtained a complete orthonormal set of 
functions {xk,r,d{y)} over the magnetic torus; those functions are labeled by a set of good 
quantum numbers (k, r, d). It was rather easy to diagonalize the magnetic Laplacian in the 
fc-space representation. Applying a unitary transformation by {xk,r,d(y)} to the eigenstate of 
the magnetic Laplacian, we finally obtained the eigenfunction in the y-spa.ce representation. 
The eigenvalues of the magnetic Laplacian were naturally interpreted as sums of energies 
of cyclic motions in the transverse directions to the magnetic field and energies of linear 
motions in the longitudinal direction to it. 

New results of this paper are the definition and representations of the magnetic algebra, 
the proof of irreducibility of the space of twisted periodic functions as a representation 
space of the magnetic algebra, the complete orthogonal set of functions ()6.8|) which provides 
a basis of the extended Fourier analysis, the eigenfunctions (17.11)1 of the magnetic Laplacian 
in explicit forms, and the eigenvalues (|7.12j) . 

Before closing this paper we would like to discuss briefly possible directions for further 
development. We treated only the Laplace operator in this paper but for application to 
physics it is more desirable to treat the Schrodinger operator 

H = ~A + V, (8.1) 

which has a potential energy term V. The potential V is a periodic function; in the t- 
coordinate it satisfies Vit 1 , . . . , f + 1, . . . , t n ) = Vit 1 , . . . , f , . . . , t n ) for each i. It acts on the 
twisted periodic function f(t) by multiplication. To take the potential term into account 
we may introduce new operators X k by 

(X*f)(t) = e 2mtJ f(t) (j = 1,2,..., 2m), (8.2) 

which belong to the same family 
operator can be expanded as 

oo 

v(t\...,n = e 

CTl,<T2,...,cr„ = 
OO 

E 

<ri,<T2,...,0n= 



of operators T of ()5.9jl . Then any periodic potential 



2wi(a 1 t 1 + -+a n t n ) 



C (A 1 ) 01 ■ ■ ■ (J^ 2m )°" 2m (7 nl )°" 2m + 1 . . . (rpn—2m\cFn 3^ 
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We can easily calculate commutators of X's with other operators to get an algebra which 
is an extension of the magnetic algebra. The resulted algebra is isomorphic to the so-called 
noncommutative torus [21] although we do not yet examine these relation thoroughly. 

Another direction for future development is to solve an eigenvalue problem of the Dirac 
operator in the n-torus in the background magnetic field. We also construct supersymmetric 
field theory, which have both scalar and spinor fields as its constituents to pursue a new 
mechanism of supersymmetry breaking. 
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Appendix A. Distribution of zeros 

As announced at (j4.11j) we prove existence of a transformation matrix L that has zeros in 
the pattern 

1 L 2 2 P \ L\ 1 i 2 2 m+ r \ / * > 
L= L\_ x L% L 2 i m+r = * * (A.l) 

j2m+k j2m+k j2m+k J \ . . . 
\ L 2p-l L 2q h 2m+r I \ * * * I 

with i,j,p,q = 1, . . . , m and k, r — 1, . . . , n — 2m. 

In R n we have an anti-symmetric bilinear form B. We say that a vector u is longitudinal 
with respect to B if it satisfies 

B(u,v) = (A.2) 

for arbitrary v G R n . We call 

M° = {ue R n | Vv G R n , B(u, v) = 0} (A.3) 

a longitudinal vector subspace. Let (t\t 2 , . . . ,t n ) be the coordinate system that expresses 
B in the standard form 

in m 

b = - <t>a dti A dtj = £ dt2j ^ A df2j - ( A - 4 ) 

as in (|3.14j) . Let (y 1 , y 2 , . . . , y n ) be another coordinate system that is related to (t 1 , t 2 , . . . , t n ) 
by a linear transformation y l = Y^=i L l j P . The matrix L is not yet specified. Basis vectors 
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generated by these coordinates are related as 



d n Bv l d n d 

dP {ri dP dy* ^ j dy v 1 ; 

d n dP d n d 

W = ^ 1 Wd¥ = ^i {L ~ 1)) dp- (A ' 6) 

Define vector subspaces M~ and M + of R n as 

M ' = R W® R W®"'® R dJ^> (A ' 7) 

M+ = R w*® R w®---® R ^- (A ' 8) 

Then M~ © M + © M° = R n . Now let us remember that R n is equipped with inner product 
structure. Then the two-form B : R n x R n — > R can be regarded as an antisymmetric 
operator B : R n — > R n . As a square of a linear operator B 2 is well-defined and becomes 
a symmetric operator and therefore is diagonalizable by an orthogonal transformation. B 2 
has non-positive eigenvalues. The eigenspace W° associated with the zero eigenvalue of B 2 
coincides with M°. Of course, 



d d d 



\ dt 2m+l ' dt 2m + 2 ' " ' ' dt n 

is a basis of M° = W°. We take an orthonormal basis 

{ d d d 



(A.9) 



(A.10) 



1 Qy2m+1 ' Qy2m+2 ' ' ' ' ' J 

of W . This implies that 

^ta+r = L im+r = (L ^^m+r = 1 ) 2m+r = (A. 11) 

in (TPJ) and (jOJ). 

Let us step into a difficult part of the proof. Each eigenspace associated with each 
nonzero eigenvalue of B 2 can be decomposed into two-dimensional subspaces such that each 
two-dimensional subspace Wj is irreducible with respect to the action of B. Thus we get an 
orthogonal decomposition 

R n = Wi JL W 2 JL ■ ■ • JL W m JL (A.12) 
Next we define vector subspaces 

iy+ = w^- n (M + © M°), (A. 13) 

WJ = Wj n (M + © M )" 1 , (A. 14) 
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then we can show that both Wt and Wj have one dimension. First, note that dimWj- + 
dimH^ - = dimHj = 2. If dimW^ = 2, Wj coincides with Wj itself. Since the two-form 
B is degenerated on (M+ © M°), it must be degenerated also on Wj = W+ C (M+ © M°). 
This contradicts the fact that Wj is irreducible with respect to B. On the other hand, if 
dim W~ = 2, W~ coincides with Wj itself. Then we can take an arbitrary one- dimensional 
subspace M_ x C Wj = W~ C (M + © M°) L . Since W,- is irreducible with respect to B, B 
is degenerated on M_i © M + © M°. This contradicts the fact that M + © M° is a maximal 
degenerated subspace of B. Hence we conclude that dimVFy" = dimH^ - = 1. 

We take a normalized vector d/dy 2 ^ of Wj~. And we take another normalized vector 
d/dy 2 ^ 1 of Wf such that Vj = B(d/dy 2j ~ l , d/dy 2j ) > for each j = l,...,m. Then 
we obtain a complete orthonormal basis {d/dy l \i = 1,2, ... ,n} that expresses B in the 
standard form B = Ylf=i v j dy 2 ^ 1 A dy 2j . 

Since d/dy 2q <G W+ C M+ © M°, we can say that 

= (A.15) 
in (|A.6[) . By an elementary argument of linear algebra we can say that 

L*- 1 = 0. (A.16) 

The proof is over. 



Appendix B. Fourier transformation of the Hermite 

polynomials 

In our convention the Hermite polynomial is defined as 

H n (Z) = (-I)" — e< (B.l) 
The formula (j7.10|) can be deduced by a partial integration and a change of variables as 

dke 2mkz ■e^ 2/v e- 27rfe > 



f 00 dke- 2 ^( 

j — oo y 



dk 



dk 



n \ n 

a 



3 71 " l dk e~ 2nk ' v ( — j e ^ k + ivz ) l v 



e 



oo 
oo 



K "-~ ' dke~ 27Tk ' 2/v f — - — ) e 7r ( fc + il/2 ) 2 /^ 
o \ v dz ) 

V 3z J J-oo 

e W23 (--^) v^e~ W . (B.2) 
\ v az I 
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